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ABSTRACT 

Using the new Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope 
(HST), we obtained moderate-resolution, high signal/noise ultraviolet spectra of 
HD 209458 and its exoplanet HD 209458b during transit, both orbital quadra- 
tures, and secondary eclipse. We compare transit spectra with spectra obtained 
at non-transit phases to identify spectral features due to the exoplanet's expand- 
ing atmosphere. We find that the mean flux decreased by 7.8 ± 1.3% for the C II 
1334.5323 A and 1335.6854 A lines and by 8.2 ± 1.4% for the Si III 1206.500 A 
line during transit compared to non-transit times in the velocity interval -50 to 
+50 km s"^. Comparison of the C II and Si III line depths and transit /non- 
transit line ratios shows deeper absorption features near -10 and +15 km s~^ 
and less certain features near -40 and +30-70 km s~^, but future observations 
are needed to verify this first detection of velocity structure in the expanding 
atmosphere of an exoplanet. Our results for the C II lines and the non-detectio n 



of Si IV 1394.76 A absorption are in agreement with lVidal-Madjar et al.l ( 120041 ). 
but we find absorption during transit in the Si III line contrary to the earlier 
result. The 8 ± 1% obscuration of the star during transit is far larger than the 
1.5% obscuration by the exoplanet's disk. Absorption during transit at velocities 
between -50 and +50 km s^^ in the C II and Si III lines requires high-velocity 
ion absorbers. Assuming hydrodynamic model values for the gas temperature 
and outflow velocity at the limb of the outflow as seen in the C II lines, we find 
mass-loss rates in the range (8-40) xlO^'' g s~^. These rates assume that the 
carbon abundance is solar, which is not the case for the giant planets in the solar 
system. Our mass-loss rate estimate is consistent with theoretical hydrodynamic 
models that include metals in the outflowing gas. 

Subject headings: planets and satellites: atmospheres — planets and satellites: 
individual (HD 209458b) — planets and satellites: physical evolution — stars: 
individual (HD 209458) — ultraviolet: stars 



"'^Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the Data 
Archive at the Space Telescope Science Institute. STScI is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations are associated with 
the HST GTO program #11534. 
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INTRODUCTION 



Mazeh et al.l fl2000[ ) monitored the radial velocities of the GO V star HD 209458 to de- 
rive the orbital parameters and minimum mass of its transiting planet HD 209458b. Very 
shortly thereafter, accurate measurements of the mass and radius of HD 209458b wer e ex- 
tracted from radial velocit y and photometric transit observations by iHenry et al.l ( 120001 ) and 
Charbonneau et al.l (120001). Table 1 l i sts the presently accepted properties of HD 209458 and 
HD 209458b cited bv iKnutson et al.l J^OoJ. HD 209458b is likely the best-studied Jupiter- 
like exoplanet located very close to its host star. As a result of its proximity, HD 209458b 
receives very strong incident radiation (see estimate in Table 1) and stellar wind flux from 
its host star. 

HD 209458b is the first t ransiting planet for which atmospheric absorption was observed 
in the resonance lines of Na I (jCharbonneau et al.ll2002l ) and subsequently in lines of H I, C II, 
O I, and others. Using the G140M grating of Space Telescope Imaging S pectrograph (STIS) 



on the Hubble Space Telescope HST with its resolution of ~ 30 km s~^, IVidal-Madjar et al. 
(120031 ) found that the Lyman-a line flux was reduced by 15 ± 4% during transit at velocities 
between -130 and +100 km s~^. Based on this first detection of H I absorption at or 
above the planet's Roche lobe, they concluded that hydrogen is escaping from the planet. 
Vidal-Madjar et al.l ( 120041 ) then used the low-resolution G140L grating of STIS to detect an 
absorption depth at midtransit of 5 ± 2% for the unresolved Lyman-a line, 13 ± 4.5% for 
the O I 1304 A multiplet, and 7.5 ± 3.5% for the C II 1335 A doublet. Their detection of 
absorption in the O I and C II lines at or above the Roche lobe confirmed significant mass 
loss from the planet by hydrodynamic outflow. They argued that reduction in the Lyman- 
a line flux at ~ 100 km s~^ from line center indicates high velocities of neutral hydrogen 
atoms. Additionally, more absorption on the blue side of the Lyman-a emission line than 
the red side indicates outflow towards the observer and away from the planet. A subsequent 
observation of HD 209458 using the HST Advanced Camera for Surveys (ACS) showed a 
reduction in the unresolved Lyman-a line flux of (8.0 ± 5.7)% during transits, consistent 
with the STIS G140L and G140M resu lts, thereby confirming mass loss from the planet's 
large exosphere ( lEhrenreich et al.ll2008l ). 



There have been a number o f reanalyse s of th e data and fu t her d is cussions of the STIS 



transit G140M observations by iBen-Jaffell ( l2007h . ISen-Jaffell ( l2008l ). IVidal-Madiar et al. 



( 120081 ) and others. There is a consensus that HD 209458b has an extended exosphere and 
is losing mass, but important details are not yet understood. Such questions as: whether 
the exosphere size is larger or smaller than the projected Roche lobe, what are the outflow 
speed and mass-loss rate, and whether or not the mass loss from the planet produces a 
cometlike tail are not yet answered with the existing STIS and ACS observations. While 
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the STIS G140M Lyman-o; spectrum obtained by IVidal-Madjar et al.l (|2003[ ) provides some 
velocity information, higher spectral resolution and signal/noise (S/N) are needed to obtain 
the velocity structure and optical depth of the outflowing gas in Lyman-a and other spectral 
lines. 

The large abundance of hydrogen and large flux in the Lyman-a line make this line a 
prime candidate for studying exoplanet outflows, but there are uncertainties (many of them 
unique) in measuring the exoplanet's mass-loss rate only from this line. These include the 
modest S/N of the existing STIS G140M grating data, possible variability of the stellar 
Lyman-a profile at the time of transit that cannot be directly measured, broad interstellar 
absorption in the core of the line that prevents measurements of low velocity absorption 
during transit, and the time- varying geocoronal and interplanetary emission in the line core 
that cannot be completely corrected for. 

Absorption of Lyman-a photons during transit over the velocity range -130 to +100 
km s~^ raises the question of the origin of the high- velocity hydrogen that is difficult to ex- 
plain by a thermal plasma at moderate temperature. iHolmstrom et al.l (120081 ) proposed that 
the high- velocity neutral hydrogen can be produced by charge exchange between stellar wind 
protons and neutral hydrogen in the planet's outflow. Such energetic neutral atoms (EN As) 
are seen in the solar system and likely occur where the stellar and planetary winds interact 
between HD 209458b and its host star. If this were the only possible explanation, then the 
Lyman-a profile during transit provides us eful information on the s tellar wind but ambiguou s 
informati on on mass loss from the planet (IHolmstrom et al.ll2008t iMurray-Clay et al.ll2009l ). 
However, iLecavelier des Etangs et al.l (l2008bl ) argued that the observed Lyman-a profile and 
the high-velocity hydrogen atoms can be simply explained by stellar radiation pressure and 
the ENA explanation required a peculiar stellar wind model. 



The Lyman-a transit observations by IVidal-Madjar et al.l (120031 . |2004| ) stimulated a 
number of interpretive papers concerning the n aechanisms for mass loss a nd lifetime of 
a Jup it er-like planet cl o se to it s host star [e.g., IVidal-Madiar et al.l (12003): iBaraffe et al. 



(|2004[ ): iKoskinen et al.l (120071): IHolmstrom et all (120081): 



drody n amic models [e.g. 
J2004h:ljaritzet al.lJ2005[) : 



lammer et al 



Koskinen et al 



(120031): ILecavelier des Etangs et all (12004 ): 



Tian et al.l ( l2005l ): lGarcia Mufiozl (l2007f ): ISchneiter et all ( l2007h : 



(120101)] a n d hy- 



Yel 



Murray-Clay et al.l ( 20091 )]. As discussed below, the theoretical models differ in their treat- 
ment of the amount of available radiative energy input (photoionization and otherwise) to 
drive the hyrdodynamic blow-off, the location of the heating, the presence of metals in the 
outflow, and the relative amount of cooling by expansion, thermal conduction, and radiation. 
These models predict very different mass-loss rates. Which of these models are consistent 
with the spectrally resolved metal lines observed during transit? 



-5- 



With the objectives of measuring an accurate mass-loss rate from the planet and testing 
different physical models for the mass loss, one needs far better data, which requires a new 
scientific instrument. One needs higher S/N spectra with higher spectral resolution and 
sufficient sensitivity to study many spectral lines formed in the extended atmosphere and 
wind of the exoplanet. In particular, it is important to use spectral lines other than Lyman-a 
to avoid unique difficulties in the analysis of this line. As we describe below, the new Cosmic 
Origins Spectrograph (COS) on HST is well designed for this task. 



2. OBSERVATIONS AND DATA REDUCTION 

Installed on the HST during Servicing Mission 4 in 2009 May, COS is a high-throughput 
ultraviolet (UV) spectrograph optimized for point sources. Descriptions of the on-orbit 
performance characteristics of COS will be presented by Green et al. (2010, in preparation) 
and Osterman et al. (2010, in preparation). During GTO Program 11534, we observed 
HD 209458 with both the G130M and G160M gratings of the COS far-UV channel to obtain 
moderate-resolution spectra covering the 1140-1790 A spectral region. In this paper, we 
analyze the spectra obtained with the G130M grating , which includes th e spectral range 
1140-1450 A. The G160M observations are presented by [France et al.l (120101 ). The COS line- 



spread function (LSF) in the far-UV is not a simple Gaussian, but it can be approximated 
by a Gaussian with a r esolution of 17,000-18,000 (16.7-17.6 km s~^) with extended wings 



(JGhavamian et al.l 120101 ). 



We have tested for the effects of the extended wings by comparing the C II 1335.6854 A 
emission line of a Cen A (G2 V) observed with the STIS E140H grating (3 km s~^ resolution) 
convolved with a Gaussian (17 km s~^ resolution) and with the COS LSF. In the line wings 
at 50-70 km s^^, the Gaussian broadens the profile by less than 1 km s^^, but the COS LSF 
broadens the profile by 2-6 km s~^. Since both the transit and non-transit line profiles are 
broadened in the same way, their ratio is not changed by the broad wings of the COS LSF. 
The broadening of the intrinsic spectral lines by the COS LSF (compared to the Gaussian 
profile) is less than 1/3 of the spectral resolution element and thus not significant. 

Table 2 summarizes the individual G130M exposures obtained at transit (near orbital 
phase 0.00), first quadrature (phase 0.25), secondary eclipse (phase 0.50), and second quadra- 
ture (phase 0.75). The total exposure times for these phases are 4946.8 s, 7941.7 s, 4956.7 
s, and 7796. s, respectively . The orbital phases listed in Table 2 are computed using the 



ephemeris of iKnutson et al.l (120071 ). Transit occurs between orbital phases 0.982 and 0.018 
( Wittenmyer et al.ll2005l ) as indicated by the optical light curve. Since the duration of transit 



is about 3.22 hours, HST covXd observe HD 209458 during the same transit on two sequential 
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spacecraft orbits. 

The S/N in the COS spectra can be lower than predicted by photon statistics as a result 
of two instrumental effects. The hex pattern formed at the intersection of the microchannel 
plate pores imposes a ±4% pattern on the spectral signal, and shadows of the wire grid 
on the detector decrease the signal as much as 20% at certain locations on the detector 
faceplate. We checked the locations of the known wire grid shadows and found that the 
spectra of the C II, Si III, and Si IV lines studied in this paper are not affected by the 
wire grids. Both effects, which are fixed at the detector, were mitigated by moving the 
grating and the resultant spectrum to different positions on the detector. In this program, 
we observed at four different grating offset positions identified by the beginning and ending 
wavelengths of the gap between the two detector faceplate elements listed in the last column 
of Table 2. In principle, division of a transit spectrum by an non-transit spectrum at the 
same grating position should cancel the signal errors introduced by both instrumental effects 
to show the true spectral difference produced by the planet in front of the star. In practice, 
the cancellation will not be perfect because of possible time variations in the stellar spectrum 
and grating mechanisms that do not return to the same position every time. The grating 
mechanism positioning errors are not fully understood at this time. We found that the 
wavelength solutions for spectra observed at different grating offset positions are not identical 
and could be off by 1-2 pixels. We registered each observed spectrum with others at the 
same grating setting by cross-correlating the interstellar feature in the C II 1334.5323 A 
line. These registrations typically involved displacements by only 1-2 pixels, far less than 
a resolution element with no significant degradation in the spectral resolution. We then 
co-added the summed spectra for the four grating offset positions using the same cross- 
correlation technique. Also, within 500 pixels of the end of each detector segment (about 
5 A for G130M), there are additional instrumental artifacts in the spectrum. For this reason, 
we did not analyze portions of the spectrum near the ends of a detector faceplate. The 
data were processed with the COS calibration pipeline, CALCOqj v2.11b (2009-09-08) and 
combined with a custom IDL coaddition procedure. 

Figure 1 shows a portion of the G130M spectra obtained during transit and all non- 
transit phases (secondary eclipse and both quadratures). The spectra are very similar, but 
there are subtle differences that provide information on the exoplanet's atmosphere and 
mass-loss rate. In this paper, we compare transit and non-transit spectra of lines of C II, 
Si III, and Si IV. The H I Lyman-a line cannot be analyzed because geocoronal emission 
through the large aperture of COS completely dominates any stellar componant. The O I 



^We refer the reader to the cycle 18 COS Instrument Handbook for more details: 



http://www.stsci.edu/hst/cos/documents/handbooks/current/cosjcover.html[ 
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lines at 1302, 1304, and 1306 A were also not analyzed because time-varying geocoronal 
emission in these lines renders comparison of transit and non-transit spectra uncertain and 
because the lines were located close to the end of a detector faceplate. 



3. RESULTS 

3.1. Si IV Line 

A critical question is whether the emission line fluxes of the host star were significantly 
different at the times when the transit and non-transit spectra were obtained, since late-type 
stars hke the Sun ofte n show time- variable emission in lines formed in their outer atmospheres 
(e.g., 



Rottmanll2006l ). Figure 2 shows the co-added transit and non-transit spectra for the 
Si IV resonance line (1393.76 A). Since the Si IV line is formed in highly ionized gas which is 
not likely to be in the atmosphere or extended wind of the exoplanet, the very similar Si IV 
fluxes during times of transit and non-transit indictate that the average flux in the stellar 
emission lines during our observations was nearly constant. The difference spectrum panel 
in Figure 2 confirms that the transit and non-transit spectra of the Si IV line are identical 
to within the noise level. 

To make the comparison of the transit and non-transit Si IV line profiles more quantita- 
tive, we compute the mean transit /non-transit flux ratio over the velocity range -50 to +50 
km s~^. Beyond this velocity range, the fluxes are less than 15% of the peak flux, and the 
S/N in each spectral resolution element is small. For the 100 km s~^ wide velocity interval 
centered on km s~^, the mean flux ratio is 0.998. To estimate the error in the mean flux 
ratio, we have computed flux ratios for many velocity intervals contained within the -50 
to +50 km s~^ velocity range that are between 60 and 95 km s^^ wide. The dispersion of 
these flux ratios indicates that the standard deviation in the mean flux ratio for the -50 to 
+50 km s~^ velocity range is ±0.014. We therefore conclude that the mean stellar fluxes 
in the other emission lines were the same during transit and non-transit times, and do not 
rescale the transit and non-transit fluxes. At positive velocities, there are regions where 
the transit /non-transit ratio may differ from unity at or above one standard deviation. We 
discuss this in Section 3.4. 



3.2. C II Lines 

The C II lines at 1334.5323 A and 1335.6854 A are bright emission lines formed in the 
chromospheres of solar-type stars like HD 209458. On October 2, 2009, we observed the C II 



lines during transit with the G130M grating at four grating positions specified by the gap 
wavelengths listed in Table 2. For each of these grating positions, there are corresponding 
observations at both quadratures and a secondary eclipse obtained at different times. 

Figure 3 shows the co-added spectra of the C II 1334.5323 A resonance line and the 
C II 1335.6854 A line. The spectra were smoothed with a 5-pixel boxcar, which is somewhat 
smaller than the on-orbit G130M spectral resolution element of 7-8 pixels. The transit and 
non-transit spectra are plotted separately and ratioed. The deep minimum observed in the 
C II 1334.5323 A line at -6.60 km s~^ ( heliocentric) is interstella r absorption due to the Eri 



cloud located within 3.5 pc of the Sun (JRedfield &: LinskyI |2008| ) . Since the radial velocity 



of HD 209458 is -14.8 km s ^, the -6.6 km s ^ interstellar absorption is centered to the red 
of the centroid of the stellar emission line. 

Both C II lines show less fiux during transit than non-transit times, although there are 
small velocity intervals where the transit /non-transit ratio exceeds unity. To test whether 
the ratios in these velocity intervals are real or noise, we have co-added the profiles of the two 
C II lines (including the portion of the 1334 A line with interstellar absorption) to get transit 
and non-transit profiles with twice the signal and higher S/N. The co-added C II profiles 
are included in Figure 3, together with difference and ratio profiles for each C II line and 
the co-added line. The co-added C II line shows absorption during transit at all velocities 
between -50 and +50 km s~^, except at zero velocity. The deep absorption features for the 
co-added C II line are interesting and will be discussed in Section 3.4. 

Using the same method that we used for the Si IV line, we find that the mean fiux ratio 
in the -50 to +50 km s'^ velocity interval is 0.924 ± 0.022 for the C II 1334.53 A line and 
0.921 ± 0.015 for the C II 1335.69 A line. For the co-added C II line, the mean fiux ratio is 
0.922 + 0.013. 



3.3. Si III Line 

Figure 2 shows the co-added transit and non-transit spectra for the Si III resonance line 
at 1206.500 A. The difference and ratio spectra show absorption at all velocities between -60 
and +60 km s~^. There are also deep absorption features which we will discuss in Section 
3.4. 

Using the same method that we used in analyzing the Si IV and C II lines, we find 
that the mean fiux ratio in the -50 to +50 km s~^ band is 0.918 + 0.014, the same within 
errors as previously seen for the C II lines. Our result for the mean absorption depth during 
transit of 8.2 + 1.4% is very different from the non-detection absorption depth of O.OI^q'q^ 
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obtained by IVidal-Madiar et al.l (|2004J) from the STIS G140L spectrum and confirmed by 
Ben-Jaffel fc Sona Hosseinil (120 lOi ). We suspect that the difference between our detection of 
Si III absorption during transit and the previous non-detection is due to time variabihty. 
Charge transfer between protons and Si"*" ions leads to a s ubstantially higher a bund ance of 



Si"'""'" than is the case for coUisional ionization equilibrium. iBaliunas fc Butler! ( 1l980l ) found 
that Si"'""'" is the dominant ionization stage in coronal plasmas at temperatures above 20,000 K 
and is important at 15,000 K. Thus, we are not surprized to see a substantial amount of 
Si"'""'" at the limb of the exoplanet's outflow where protons from the solar wind and planetary 
escape are present to ionize Si"'", but the amount of Si"'""'" could vary appreciably with changes 
in the stellar wind, planetary mass-loss rate, and temperature of the gas escaping from the 
planet's exosphere. 



3.4. Is There Velocity Structure in the Planet's Mass Loss? 

Until now we have discussed the mean transit/non-transit flux ratios for the entire line 
cores between -50 and +50 km s~^, but there are indications of departures from these mean 
ratios which indicate that the absorption in the planet's mass loss depends on velocity. To 
test this idea, we compare in Figure 4 the difference and ratio profiles for the Si III and 
co-added C II lines. We find deep absorption features near -10 and +15 km s^^ in both 
lines that are larger than the indicated errors shown in Figures 2 and 3. Since the C II and 
Si III lines are likely formed in the same region of the exosphere, we include in Figure 4 
the co-addition of the two profiles with typical error bars for comparison with theoretical 
models. The increase in S/N with the co-addition of the C II and Si III spectra allows us to 
make more definitive statements about the velocity structure. The dips near -10 and +15 
km s~^ are more than twice as large as the errors per spectral resolution element (about 17 
km s~^) and are most likely real. In addition, the ratios near -40 and +30-70 km s~^ are 
also low, but the line fluxes decrease rapidly away from line center and these low ratios are 
less certain. Future observations with higher S/N are needed to confirm this first detection 
of velocity structure in the expanding atmosphere of an exoplanet. We also note that there 
are absorption features near +20 and +40 km s~^ in the transit /non-transit ratio for the 
Si IV line which are greater than the noise level. Future observations are also needed to test 
the reality of these features. 
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3.5. Does HD 209458b have a Comet-like Tail? 



Schneider et al.l (Il998[ ) proposed that giant exoplanets located close to their host star 



will have comet-like tails of ionized gas produced by the interaction of the stellar wind and 
ions excaping from the planet's exo sphere, and that th e se tai ls could produce absorption 
when observed in front of the star. IVidal-Madjar et al.l ( 120031 ) then predicted that stellar 
Lyman-a radiation pressure on hydrogen atoms escaping from the planet's exosphere would 
form a comet-like tail trailing the planet. The presence of gas in an exte nded tail should 
produ ce more absorption at and beyond egress than during ingress phases. ISchneiter et al. 
(120071 ) calculated possible cometary wakes and the time evolution of Lyman-a absorption 
from ingress to beyond egress. We searched our C II and Si III spectra for line profile 
differences between the ingress and egress phases to see whether these ions could be present 
in significant amounts in a tail, but we found no differences larger than the noise level. Future 
observations with higher S/N are required to verify the predictions of comet-like tails. 



MASS-LOSS RATE IN A THERMAL PLASMA OUTFLOW 



The ^ 8 ± 1% obscuration of the star in the C II and Si III lines during transit far 
exceeds the 1.5% obscuration by the exoplanet's disk, indicating absorption by an extended 
atmosphere or mass loss that is optically thick in the C II and Si II I lines. Our resu lt 



for the C II hues is consistent with the 7.5tl'l% absorption depth that IVidal-Madjar et al. 



-3.4' 



( 120041 ) obtained for the unresolved C II lines, but very different from the absorption depth 



of 0.0~^a'n% that they obtained for the Si III line. Their absorption depth of 0.0]'^n'n% for the 



-0.0' 



-0.0' 



Si IV line is consistent with our result. 



Since the size of the obscuri ng material (^optically thick in thes e lines) is similar to the size 
of the exoplanet's Roche lobe (JBen-Jaffel fc Sona Hosseinill2010l ). one possible explanation 
for the transit data is that the Roche lobe is filled by mass loss from the planet. Another 
possible explanation is that the exoplanet's escaping gas forms an extende d cometary tail 



blown out by the stellar wind, as suggested by hydrodynamical simulations ( ISchneiter et al. 
2007h . 



We estimate the mass-loss rate from the planet's atmosphere by assuming that the 
outflow forms a spherically symmetric envelope around the planet which is optically thick 
in the C II resonance line covering about 8% of the stellar surface. The spherical symmetry 
assumption is a good approximation beca use the shape of the tidally indu ced Roche lobe 
as seen during transit is close to a sphere (JLecavelier des Etangs et al.l 120041 ). The envelope 
becomes optically thin at its "limb" because gas expansion and mass flux conservation reduce 
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the gas density with increasing radial distance (r) from the planet. We define the envelope 
"limb" to be at the radial distance from the planet for which the absorption area is the same 
as that of an opaque sphere of the same size. 

We consider a line of sight from the star that passes through the envelope with p the 
point of closest approach to the planet. The quantity x measures the distance along this 
line of sight from point p back to the star. The optical depth to the observer in the C II 
resonance line along this line of sight at the apparent limb of the envelope is, 

/■oo 

r(p, Az/) = 2k nc n{x)dx = 0.69, (1) 

Jo 

where nc n is the number density of C+ ions, k = -^^^^-^f- is the line opacity, / = 0.128 
is the oscillator strength, m is the electron mass, and Aud is the Doppler width. The value 
of p depends somewhat on the displacement from line center, Au, because for expanding 
gas, most of the absorption occurs near the frequency corresponding to the line-of-sight 
component of the expansion velocity. The limb occurs at r = 0.69, where as much stellar 
fiux is transmitted through the envelope for lines of sight passing inside of p as for those 
lines of sight passing outside of p, JJ^g e~^dt = J^ ' e~^dt, where t is the optical depth 
through the envelope for lines of sight with different values of p. Since the gas density and 
thus optical depth increase along lines of sight inside of p, the limb will occur at somewhat 
smaller optical depths depe n ding o n details the density and outfiow velocity. For example. 



Lecavelier des Etangs et al.l (l2008af ) calculated r = 0.56 for an atmosphere in hydrostatic 
equilibrium, which has a density that decreases exponentially outward but does not have an 
outflow or heating in the outer layers. If the density decreases with radial position according 
to a power law, then r could be smaller than 0.69 by a factor of 2 or 3. The resulting mass 
loss rates would also be smaller by this factor. We call attention to this uncertainty that 
should be with resolved by realistic models based on physical principles that predict transit 
line profiles consistent with the COS data. 

Since the gas density decreases with increasing radial distance from the planet r (thus 
increasing x along the line of sight), p should decrease slightly with incre asing Au. We 



compute Az/n = 2.7 x 10 " s , assuming a gas temperature of 10,000 K (I Garcia Munoz 



20071 : iBallester et al.ll2007l ). If the outfiow is turbulent or some absorbers have superthermal 
velocities (see below), then Aud will be larger. For a spherical outfiow with constant mass 
fiux, the mass fiux of C^ ions from the planet is 

Mc II = A'nr'^mcvnc ii{r), (2) 
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where v is the outflow speed at r, mc is the mass of a carbon atom, and r^ = x^ + p^. 

We estimate Mq h by requiring that the gas in the exosphere obscure 8% of the stellar 
disk within Az//? of line center, corresponding to ±3.7 km s^^. This condition will occur 
when the optical depth is about 0.69 at p^ = 0.08i?^, i.e., 



r(j9 = 0.284i?^, Az/) 



2^e''fMc II 



The integral has an analytical solution by setting u = x/p of 



'^ dx 

-r = 0.69. 



(3) 



dx 



du 



IT 



P Jo M + 1 2p 



(4) 



We adopt t | = 1 km s ^ at p = 0.28 4-R^ = 2.36Rpiane.t from the th eoretical models of 
Tian et al.l feoosl ). ICarcia Munozl fcoOTJ ). and iMurrav-Clav et al.l fl2009h . We then find that 



Mc II = 2.1 X lO^g s"^ = 6.8 X lO^^g yr'^ 



(5) 



In his model for the exosphere of HD 209458b. iGarcia Munoa (120071 ) found that essentially all 
of the carbon is sii igly ionized. If we as sume that the carbon abundance is solar, 2.7 x 10~^ 
that of hydrogen (jAsplund et al.ll2009l ). the total mass- loss rate from the planet for only 
thermal broadening of the C II lines is 



M 



total 



-1 



18, 



-1 



8.0 X lO'^g s"' = 2.4 X 10'«g jT-\ 



(6) 



Since the planet's mass is 1.2 x 10^° g (JKnutson et al.l 120071 ) . the fractional mass-loss 
rate per year is very small, as shown by 



Mtotal 



M 



2.0 X lO^^V^^- 



(7) 



planet 



This calculation assumes that the mass-loss rate and planetary orbit do not change with 
time. Since additional line broadening is required to explain the broad absorption in the 
C II lines, these values of Men and Mtotai are at the lower end of their realistic ranges (see 
Section 5). 
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DISCUSSION 



We now compare our empirical mass-loss rate with those of various theoretical models. 
All theoretical models of which we are aware assume hydrodynamical outflow driven by high 
temperatures and pressures near the base of the outflow. The important questions are the 
amount and location of energy absorbed from the star, the important cooling mechanisms, 
and the location of the sonic point in these transonic (Pa r ker-lik e) winds relative to the 
location of the exoplanet's Roche lobe. IVidal-Madjar et al.l (120031 ) demonstrated that neu- 
tral hydrogen must be escaping from a volume larger than the planet's Roche lobe by an 
unspecified mass-loss mechanism and that stellar radiation pressure further accelerates the 
outflow and shapes a comet-like tail. They proposed that the mass-loss rate is M > 10^° 
g s~\ but they added that "owing to saturation effects in the [Lyman-a] absorption line, a 
flux larger by seve r al ord ers of magnitude would produce a similar absorption signature." 
Murray-Clay et al.l (120091 ) computed models with and without radiation pressure and con- 
cluded that radiation pressure does not increase M, but it can change the location of the 
mass loss from the day side to the night side of the exoplanet. Tidal forces also do not 
change M apprec iably, but they move the sonic point closer to the exoplanet and increase 
the flow velocity JCarcia MuiiozlboOTl : iMurrav-Clav et al.ll2009[ l. 



All of the theoretical models assume that HD 209458 is a main sequence star like the Sun 
with ultraviolet and extreme ultraviolet radiative output similar to the Sun at low or high 
activity. Since mass-loss rates are proportional to the energy input (the so-called energy- 
limited escape rate), the critical question is the energy input rate. In all of the models, 
expansion (PdV work) is the domin ant cooling mechanism, although radiation by H|^ at the 
base of the flow must be included (lYellell2004l ). and radiation in the Lyman-a line can be 
important for T Tauri star radiative input rates. Many models assume that only the energy 
remaining after stellar Lyman continuum (E> 13.6 eV) photons ionize neutral hydrogen is 
available to energize the outflow, so that the assumed fraction of this radiation (83% or 100% 

For 
4.7 X 10^° ' ^^ "" 



in differe nt models) c ontrols the outfl ow. For such models, M is a few ti mes 10^° g s ^ 



example, lYelld (120041), as corrected in JYellel (120061 ). finds M 



find M < 6 X 10^° 

l ow ac tivity Sun but M 

J2009h find that M 



Tian et all (I2005D 



Garcia Munod (120071 ) finds M = 6.1 x 10^" for EUV flux similar to the 



1.5 X 10^^ for EUV flux like the active Sun. iMurray-Clay et al. 



2 X 10^° g s-^ 



The inclusion of atoms and molecules other than H, H"*", and H2 can significantly change 
the outflow, in particular by increasing the mass-loss rate, because more of the stellar UV 
radiation is available to heat the outflow and drive the mass l oss. For example, the Lyman- 
a line contains more flux than the entire Lyman continuum (JRibas et al.ll2005l ). and other 
UV emission lines and continuum radiation in the wavelength range 912-2000 A can be 
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absorbed i n the bound-fr e e con tinua of atomic species (e.g., CNO) and molecules including 
H2. When I Garcia Munozl (120071 ) added solar abundances of CNO to his models, M increased 
by a factor of 4.5 (see his Table 6). Including partial absorption of the 1000-2000 A stellar 
flux by atoms and molecules with solar abundances increases M to 6.3 x 10^^ g s~^. Metals 
must exist in the distant outflow since we see absorption during transit in the C II and Si III 
lines. By equating the energy input of both Lyman c ontinuum and Lyman-g with stellar 
fluxes estimated from the rotation rate of HD 209458, iLecayelier des Etangsi (120071) found 
M ~ 3 X 10^^ g s^^. Usin g a similar e nergy balance argument, iLammer et al.l ( 20031 ) obtained 



M ~ 3 X 10 g s , but lYelld (l2004j ) suggested that this result could be an overestimate as 



a result of not including Hjj' cooling and ionization of H. 

A critical question is whether line absorption during transit is produced entirely by 
the outflowing thermal plasma or whether high-velocity atoms and ions (either at higher 
temperatures or nonthermal) are required to e xplain the broadest part of the absorption 
line profiles. iBen-Jaffel fc Sona Hosseinil ( I2OIOI ) show that neutral hydrogen absorption by 
scaled hydrodynamic outflow models may explain the flux decrease during transit even at 
±100 km s~^ from the center of the Lyman-a line, but such models predict much less 
absorp tion during transit than is detected in the unresolved O I and C II line flux measure- 
ments (jVidal-Madjar et al.l 120041 ). Additional high- velocity absorbers are needed to explain 
the O I, C II, and probably also the Lyman-a data, but there is no agreement concern- 
ing the source of the hig h-velocity absorbers. One possibility is radiation pressure, but 
Murray-Clay et al.l ( 120091 ) argue that the observed absorption near +100 km s~^ from the 
center of the Lyman-g line is difficu lt to explain by stellar radiation pressure as originally pro- 
posed bv lVidal-Madiar et al.l (120031). However, radiation pressure on neutral hydrogen atoms 
must be present, and lLecavelier des Etangs et al.l (l2008br) argue that radi ation pressure by it 



self ca n explain the Lyman-a observations. iHolmstrom et al.l (120081 ) and iMurray-Clay et al. 



(120091 ) argue that the observed decrease in Lyman-a flux near and beyond ±100 km s~^ 
during transit requires high-velocity H atoms that could be explained by charge exchange 
between stellar wind protons and neutral H atoms in the exoplanet's outflow. These high- 
velocity H atoms are called EN As (energetic neutral atoms). 



Ben-Jaffel &: Sona Hosseinil (120101 ) have tested whether the presence of high-velocity 



atoms and ions could explain the flux reduction observed in the unresolved O I and C II 
line transit data. They assume that a portion of the gas has an effective temperature 
(either thermal or nonthermal) many times larger than the thermal temperature in outflow 
models. For the unresolved C II 1334, 1335 A multiplet, they can explain the observed flux 
reduction by C"^ ions with effective temperatures 5.5-107 times larger than the background 
gas temperature {Tb) in typical hydrodynamic models depending on the properties of the 
outflow models and the assumed location of the C"*" ions. In Figure 9 of their paper, they 
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show examples of predicted C II 1334 A line profiles during transit with Tcii/Tb ~ 16 and 
~ 28 that can explain the total decrease in unresolved line fiux during transit, but these 
particular examples are either too narrow or triangular in shape, unlike the observed C II 
line profiles which show significant flux reduction between -50 and +50 km s~^ with probable 
enhanced absorption near -10 and +15 km s~^. Ben-Jaffel (private communication) informs 
us that other models with somewhat larger Tcii/Tb ratios can explain the flat reduction 
in C II line flux between -50 and +50 km s~^. However these models have ad hoc velocity 
distributions and they do not explain the observed velocity structure. In their paper, they 
conclude that COS spectra are needed to "...reveal the true balance between thermal and 
non-thermal populations..." . We encourage the development of new models, preferably based 
on plausible physical mechanisms, to explain the observed line proflles. 

Our analysis of the C II line absorption during transit assuming only thermal absorption 
by 10,000 K gas and solar carbon abundance led to our estimate of Mtotai = 8.0 x 10^° g s^^. 
This mass-loss rate is somewhat larger than the theoretical estimates mentioned above that 
assume an outflow with no metals. The assumption of no metals in the outflow is clearly 
unrealistic as we observe absorption by C II and Si III at 2.36Rp/anet- At the same time, our 
analysis of the C II lines assuming only thermal line broadening (at 10,000 K) is unrealistic 
as the thermal velocity is only 3.7 km s~^, and we observe absorption out to ±50 km s~^. 
Since we measure the mass flux at the limb of the outflow where the C II optical depth 
is small (r = 0.69), we cannot invoke optical thickness to explain the broad proflle. Note 
that the same is true for th e Lyna an-g line, but (as the referee has called to our attention) 



Ben-Jaffel &: Sona Hosseinil ( 120101 ) and other authors include opacity and Lorentzian line 



shapes at the limb as a possible mechanism for absorption line broadening. 

Equation 3 shows that Men is proportional to the Doppler broadening parameter Aud. 
To obtain 8% absorption at ±50 km s~^ from line center requires either turbulent broadening 
with velocity of about 50 km s~^ or superthermal C"*" ions with similar velocities. If turbulent 
broadening is responsible, then the mass-loss rate increases by a factor of 13 to Mtotai = 
1 X 10^^ g s~^. However, 50 km s~^ turbulent velocities are highly supersonic and unlikely. 
We estimate that a sensible upper limit to the turbulent velocity would be twice the flow 
speed or 20 km s~^, for which Mtotai ~ 4 x 10"^^ g s~^. This would increase Mtotai / M^ianet to 
TnAx 10"^^ yr~^ 

If high-velocity C"*" ions are responsible for the broad absorption, then the mass-loss 
rate could be closer to the thermal value because only a small fraction of the C"*" ions are 
likely to have high velocities. We conclude that Mtotai = ( 8-40) x 10^° g s~^ is c onsistent 



with the C II line proflles. We flnd it interesting that when I Garcia Munozl (120071 ) includes 



metals with solar abundances in his hydrodynamic mass outflow models, the mass-loss rate 



-16- 



increases from M = 1.4 x 10^^ to M = 4.95 x 10^^ g s~^, a value near the top end of our 
empirical mass-loss rate estimate. Thus high- velocity C"^ ions are required to explain the C II 
absorption during transit. Similar arguments would likely explain the broad Si III transit 
absorption. 



6. CONCLUSIONS 

The high sensitivity and moderate spectral resolution of COS allowed us to obtain the 
first measurements of gas absorption in the exosphere of a transiting exoplanet in lines of 
C II and Si III with velocity resolution. During transits of HD 209458, we found that the 
mean flux in these lines is reduced by 7.8 ± 1.3% for the co-added C II lines and 8.2 ± 1.4% 
for the Si III line in the velocity range -50 to +50 km s~^. The 8% absorption constraint 
sets the location of the outflow limb {2.36Kpianet) where the C II lines become optically thin. 
Theoretical hydrodynamic models show that the outflow velocity is about 10 km s^^ and 
gas temperature about 10,000 K at this location. If we assume that the outflow is purely 
thermal, then Mtotai = 8 x 10^° g s~^, somewhat larger than theoretical estimates based on 
outflows with no metals. Since absorption in the C II lines is much broader than indicated 
by Doppler broadening with a thermal velocity of 3.7 km s~^, the lines must be broadened by 
large turbulence or high-velocity C"*" ions. If turbulent broadening dominates, then the 8% 
absorption seen out to ±50 km s~^ would require mass-loss rates as large as 1 x 10^^ g s~^, 
but such highly supersonic velocities are not likely. Instead, absorption by high-velocity C"*" 
and Si"'""'" ions formed by the interaction of the stellar wind and planetary mass loss or some 
other mechanism broadens the C II and Si III lines. We conclude that the mass-loss rate 
in the exoplanet's outflow is in the range (8-40) x 10^° g s^^. These rates assume that the 
carbon abundance is solar, which is not the case for the giant planets in the solar system. 
We note that mass-loss rates predicted by hydrodynamic models that include solar abundant 
metals in the outflow lie at the upper limit of this range. 

Comparison of the C II and Si III line depths and transit/non-transit line ratios shows 
deep absorption features near -10 and +15 km s~^ in both lines. Since the C II and Si III 
lines are likely formed in the same region of the exosphere, we include in Figure 4 the co- 
addition of the two proflles with typical error bars for comparison with theoretical models. 
The increase in S/N with the co-addition of the C II and Si III spectra allows us to make 
more deflnitive statements about the velocity structure. The dips near -10 and +15 km s~^ 
are more than twice as large as the errors per spectral resolution element (about 17 km s^^) 
and are most likely real. In addition, the ratios near -40 and +30-70 km s~^ are also low, but 
the line fluxes decrease rapidly away from line center and these low ratios are less certain. 
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Future observations with higher S/N are needed to confirm this first detection of velocity 
structure in the expanding atmosphere of an exoplanet. 

The COS transit observations of the absorption velocity structure in these fines demon- 
strate tfie need for including fiigfi-velocity C"*" and Si^^ ions in models of tfie exospfiere 
of HD 209458b and presumably otfier close-in exoplanets. It remains to be seen wfietfier 
tfi e inclusion of a simple distr i bution of ions witfi fiigfi effective temperature as proposed 



by iBen-Jaffel fc Sona Hosseinil (J20ld ) is adequate to explain tfie data. Tfie new COS data 
provide tfie opportunity to test a range of absorption models based on sucfi pfiysical mecfia- 
nisms as radiation pressure, cfiarge excfiange, and interactions between tfie stellar wind and 
tfie planet's outflow to understand mass- loss processes in tfie exospfieres of close-in planets. 

Tfiis work is supported by NASA tfirougfi grants NNX08AC146 and NAS5-98043 to 
tfie University of Colorado at Boulder. We tfiank tfie many people at tfie Space Telescope 
Science Institute, Goddard Space Fligfit Center, Ball Aerospace and Tecfinologies Corp. and 
CASA at tfie University of Colorado for tfie excellent fiardware and software tfiat fias made 
COS a great success. We particularly wisfi to tfiank tfie referee for fiis many tfiougfitful and 
tfiorougfi comments and fiis bringing to our attention important publications tfiat fiave led 
to significant modifications in tfie paper. We tfiank Dr. J. Scfineider for encouraging us to 
look for an assymetry between ingress and egress. We also tfiank Dr. Lofti Ben-Jaffel for 
clarifying our understanding of tfie Ben-Jaffel & Hosseini (2010) paper. We tfiank Dr. Alain 
Lecavelier for calling our attention to tfie different values of r at tfie limb tfiat depend on 
tfie density structure in tfie envelope. 
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Fig. 1.— A portion of the COS G130M spectrum of HD 209458. The lower (black) plot is 
the sum of four spectra obtained during transit. The upper (red) plot (displaced upward) is 
the sum of the ten quadrature and four secondary eclipse spectra. 
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Fig. 2. — Upper panels: Comparison of the four co-added spectra obtained during transit 
(black) with the co-added non-transit spectra (red) for the Si III and Si IV lines. Repre- 
sentative error bars per pixel are included. Middle panels: difference spectra (transit minus 
non-transit) with representative error bars. Lower panels: flux ratio (transit /non-transit) 
spectra. The ratios are not plotted beyond ±50 km s^^ as the errors are large at low flux 
levels. The velocity scale is relative to the -14.8 km s~^ radial velocity of the star. 
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Fig. 3. — Upper panels: Comparison of the four co- added spectra obtained during transit 
(black) with the co-added non-transit spectra (red) for the C II 1334.53 A and 1335.69 A 
lines. The upper right panel show a co-addition of both C II lines on the same velocity 
scale. Representative error bars per pixel are included. Middle panels: difference spectra 
(transit minus non-transit) with representative error bars for each C II line and the combined 
C II line. Lower panels: flux ratio (transit /non-transit) spectra. The ratios are not plotted 
beyond ±50 km s~^ as the errors are large at low flux levels. The velocity scale is relative 
to the -14.8 km s~^ radial velocity of the star. 
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Fig. 4. — Upper and middle panels: Comparison of the difference spectra (transit minus non- 
transit) and ratio spectra (transit /non-transit) for the co-addition of the C II Lines (black) 
and the Si III line (red). Lower panel: Co-addition of the C II and Si III ratio spectra with 
typical errors per spectral resolution element (about 17 km s~^). The horizontal line at 0.92 
is the mean ratio for the velocity interval -50 to +50 km s^^. The dips near -10 and +15 
km s~^ are each more than 2a and likely real. The low ratio features near -40 and +30-70 
km s~^ occur where the line fluxes are low and are thus less certain. 
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Table 1. Properties'' of HD 209458 and HD 209458b 



Spectral type 



GOV 



Distance (pc) 


47 


M,/Mq 


-, -|Q. +0.066 


Ri,/ Rq 


1 1 OCT +0.020 
±.±ZO_QQ23 


^orbit (days) 


3.52474859(3 


-n^planet / -'*' Jup 


1 Q9n+0-024 


Mpianet/Mjup 


0.64 ±0.06 


Pplanet/Pjup 


0.26 ±0.04 


Semimajor axis (AU) 


0.045 


Stellar flux at planet/Jupiter 


13,400 


Escape speed (km/s) 


42 ±4 


Orbital Velocity ampl. (km/s) 


146 


Transit duration (lir) 


3.22 


Transit depth for -Rpianet 


1.5% 


Transit depth for -RRoche 


10% 



''Data from iKnutson et al.l (120071 ). 
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Table 2. COS G130M Observations of HD 209458 (Program 11534) 



Exp. Time (s 


) Phase 


Day 


Start 


Stop 


Gap' 


'(A) 


2340.2 


0.2516-0.2593 


19 Sep 


10:10:15 


10:49:15 


1278- 


-1290 


0955.2 


0.2684-0.2716 


19 Sep 


11:35:37 


11:51:32 


1288- 


-1299 


1851.2 


0.2722-0.2782 


19 Sep 


11:54:40 


12:25:31 


1296- 


-1306 


0560.0 


0.2873-0.2891 


19 Sep 


13:11:29 


13:20:49 


1288- 


-1299 


2235.1 


0.2898-0.2971 


19 Sep 


13:24:06 


14:01:21 


1306- 


-1316 


2340.2 


0.7059-0.7136 


24 Sep 


13:12:09 


13:51:09 


1278- 


-1290 


0945.2 


0.7229-0.7260 


24 Sep 


14:38:04 


14:53:49 


1288- 


-1299 


1787.2 


0.7266-0.7324 


24 Sep 


14:56:57 


15:26:44 


1296- 


-1306 


0925.2 


0.7417-0.7448 


24 Sep 


16:13:54 


16:29:19 


1288- 


-1299 


1798.2 


0.7454-0.7513 


24 Sep 


16:32:36 


17:02:34 


1306- 


-1316 


1045.2 


0.9913-0.9947 


02 Oct 


14:32:10 


14:49:35 


1278- 


-1290 


1096.2 


0.9954-0.9990 


02 Oct 


14:52:52 


15:11:08 


1288- 


-1299 


1400.2 


0.0077-0.0123 


02 Oct 


15:55:28 


16:18:48 


1296- 


-1306 


1405.2 


0.0129-0.0175 


02 Oct 


16:21:56 


16:45:21 


1306- 


-1316 


1057.2 


0.4880-0.4915 


18 Oct 


10:57:04 


11:14:41 


1278- 


-1290 


1093.1 


0.4921-0.4957 


18 Oct 


11:17:49 


11:36.02 


1288- 


-1299 


1401.2 


0.5043-0.5089 


18 Oct 


12:19:37 


12:42:58 


1296- 


-1306 


1405.2 


0.5095-0.5141 


18 Oct 


12:46:06 


13:09:31 


1306- 


-1316 



^The gap is the wavelength interval between the two detector faceplate el- 
ements. No data are available for the gap region, and there are instrumental 
artifacts at wavelengths close to the gap. 



